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We develop the cluster expansion and the Mayer expansion for the self-gravitating thermal gas 
and prove the existence and stability of the thermodynamic limit N, V — > oo with N/V~s fixed. 
The essential (dimensionless) variable is here r\ = Crr [ — (which is kept fixed in the thermodynamic 

1/3 T 

limit). We succeed in this way to obtain the expansion of the grand canonical partition function in 
powers of the fugacity. The corresponding cluster coefficients behave in the thermodynamic limit 
as (jS/Y Cj where Cj are pure numbers. They are expressed as integrals associated to tree cluster 
diagrams. A bilinear recurrence relation for the coefficients Cj is obtained from the mean field 
equations in the Abel's form. In this way the large j behaviour of the Cj is calculated. This large 
j behaviour provides the position of the nearest singularity which corresponds to the critical point 
(collapse) of the self-gravitating gas in the grand canonical ensemble. Finally, we discuss why other 
attempts to define a thermodynamic limit for the self-gravitating gas fail. 
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I. INTRODUCTION 

The self-gravitating gas has been the subject of attention since many years 0, 0, EJ 01 I n refs.0,0] we recently 
investigated the self-gravitating thermal gas using Monte Carlo simulations, mean field methods and low density 
expansions Q. We have shown that the system possess a well defined infinite volume limit in the grand canonical 
(GC), the canonical (C) and microcanonical (MC) ensembles when N, V —> oo keeping N/Vs fixed. A relevant 
variable here is dimensionless ratio 

Gm 2 N , . 

"=173^-. (1-1) 
V 3 2 

which is kept fixed in the N, V — > oo limit. All physical quantities per particle turn out to be functions of the single 
variable rj and are well defined and finite in the thermodynamic limit N, V — > oo with r\ fixed. 

In this paper we develop the cluster expansion and Mayer's approach to the self-gravitating thermal gas and provide 
a rigorous demonstration of the existence of the thermodynamic limit N, V — > oo with r\ fixed. 

The cluster expansion is apowerful tool allowing to express the partition function in a power series of the density 
for short range interactions [6| . We apply and adapt this method (the Mayer expansion) which is purely combinatorial 
to the self-gravitating gas. We succeed in this way to obtain the coordinate partition function Qat as a power series 



2 



in the thermodynamic limit. This is derived by generalizing the saddle point method used in ref.|6J. The expansion 
is obtained in terms of the coefficients c n which are pure numbers. More explicitly, 

i logQjvfo) N = l i g(v t n ) - log*, - l + O C^j . (1.2) 



where 



+00 



,(x) = J2cjX j - (1-3) 
j'=i 



and is the solution of the equation, 



+00 



t v g'{t v v) = 1 > i- e. ^,3 Cj (t^Tiy =r) . 

3=1 

Moreover, — -A. log(5Ar(?7) is the free energy of the self-gravitating gas minus the free energy of an ideal gas divided by 
NT. The series for 17(77 t v ) in eas. (|1.2|l - (|1.3|l is therefore a high temperature or low density expansion [see ea . (|1 . . 

The coefficients c„ can be expressed in the thermodynamic limit as a sum of 3n-uple integrals associated to tree 
cluster diagrams. Loop cluster diagrams are subdominant for N — ► 00. The coefficients c n only depend on the 
geometry of the box. The first c„ are obtained by explicit evaluation of the cluster integrals. We have for the sphere, 



_ 1 sphere 3 / 47T \ sphere ^1 f sphere 373 



Moreover, we use the connection with the mean field approach to obtain a nonlinear recurrence relation for the 
coefficients ctf here : 



c R 



1 n— 1 

l)( n _l) E^ c ^ s * (2n-2« + l) for n>2, 



(2n+l)(n-l)^ 



where 



(The label i? stands for the spherical geometry). This allows to compute systematically all and to find their large 
n behaviour as 



£ 0.309360 . 



n 2 
(ugcY 



1 + C'7 I - 



Here, wqc = 0.30034... is the radius of convergence of the series ea. ()1.3|> . Therefore, u — uqc is the nearest 
singularity of gniu) in the u-planc. 

We show below that the function g is related to the partition function Zcciili z ) m the grand canonical ensemble 
by 

N 

\ogZ GC {v,z) = — g(v h) i 

V 



where 

Out (lT\ ( m.T\ 2 



Z GC (r,,z)=^ 

N=0 



fmTY 



In addition, the fugacity {z — ct) turns out to be given by z — t n e t ( where \i and /xo are the chemical potentials 
for the self-gravitating gas and for the ideal gas, respectively. 
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The point uqc — 0.30034 . . . where g(u) is singular corresponds to the critical point (collapse) in the grand canonical 
ensemble. Recall that the collapse point depends on the ensemble considered [2| ■ 

In conclusion, our investigation here on the cluster expansion and the Mayer expansion for the self-gravitating gas 
shows that: 

• The self-gravitating gas admits a consistent thermodynamic limit N, V — » oo with — r fixed. In this limit, 

extensive thermodynamic quantities like energy, free energy, entropy are proportional to TV. That is, in this 
limit the energy, free energy, and entropy per particle are well defined and finite. 

• The cluster expansion and the mean field approach provide the same results in the thermodynamic limit 
N, V -> oo with A- fixed. 

• The partition function needs a small short-distance cutoff a in order to be well defined, the thermodynamic 
limit has a finite limit for a — ► 0. In other words, the contributions to the partition function which diverge 
for a —> are subdominant for TV — > oo (see sec. II). For large N and fixed cutoff a the potentially divergent 
contributions for a — > are suppresed by a factor at least rf 1 /N 2 compared with the dominant contribution 
for N — > oo. That is, the N, V = oo limit of the cutoff model (with fixed) has a finite limit for a — > 0. 

Therefore, subdominant corrections in r]/N can always be neglected. Realistic models of the self-gravitating 
gas (interstellar medium, galaxy distribution) require a small non-zero short distance cutoff since molecular and 
atomic forces dominate over gravitational forces for short distances. 



II. THE CLUSTER EXPANSION FOR THE SELF-GRAVITATING GAS 



We investigate in this section the self-gravitating gas in thermal equilibrium at temperature T = j3 . That is, we 
work in the canonical ensemble where the system of N particles is in contact with a thermal bath at temperature T. 
We assume the gas being on a cubic box of side L. 

The partition function of the system can be written as 



where 



N ,2 



G is Newton's gravitational constant. 

At short distances, the particle interaction for the self-gravitating gas in physical situations is not gravitational. 
Its exact nature depends on the problem under consideration (opacity limit, Van der Waals forces for molecules etc.). 
We shall just assume a repulsive short distance potential, that is, 

v A (\q l -q j \) = - - ={ (2.3) 

where A << L is the short distance cut-off. 

The integrals over the momenta pi, (1 < I < N) in ea. l|2.1|l can be computed immediately. 

It is convenient to introduce the dimensionless variables r;, 1 < I < N making explicit the volume dependence as 

qi = Lfi , n = (xi,yi, zi) , 

< xi,y h zi < 1 . (2.4) 
That is, in the new coordinates the gas is inside a cube of unit volume. 
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The partition function takes now the form, 

T\ = 



3JV 

2\ — 



1 fmTT 2 

Z C (N,T) = — {——) Q N ( V ), (2.5) 



where 



' N 



r l r l N 

{rf)= ... Y[d 3 n e W*.-.*W ; (2.6) 



o 



1=1 



is the dimensionless variable 



and u(fi, . . . , rjv) is defined by 



^(fi,...,^) = ^ E i~ . a = A/L<^l (2.8) 

l<l<j<N 1 Jl 

In this way all dependence on the volume V = L 3 is buried in the variable rj. 
The coordinate partition function Qn(v) can be written as 

/•i /■! N n _ 
Qn(v)= ■■■ Ud 3 n J] e- '^>» . (2.9) 

Jo Jo l=1 i<i<j<N 

We are now interested to expand Qn(v) m powers of In order to do this is convenient to define: 

fl . = e ^M% _ i . 

For small -3= and fixed a > we have, 



flj = ATl ^ ~, + o 



m • 

All the integrals over 77 in ca. H2.6|) are finite provided we keep a > 0. 

We can now multiply out the products of /'s in the coordinate partition function Qn(v)i 

I] (1 + fij) = 1 + E fa + E fa /« + ••■ 

i<l<j<N 

Thus, by introducing the /y functions the effects of interparticle forces are better exhibited. 

A systematic treatment of such sum of products can be found in rcf. 6] . The outcome is that a general term in the 
partition function Qn (v) ec L- 1|2.9|) , can be factorized as the product of several integrals over the coordinates fj . Each 
integral corresponds to a 'cluster' of particles and is called bj(rj, N) where 

M'/-*v.) = - / / ... tfid'n (2.i.:i) 



with 



Sj = l 



S\2 — JX2 

<Sl23 — 111 /23 + /l2 /l3 + /l3 SlZ + fl2 fl3 ^23 
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#1234 = /12 /23 hi + nineteen permutations+ 
+/12 /13 /24 /34 + fourteen permutations+ 

+/12 /13 /14 /24 /34 + five permutations + /12 /13 fu hz hi hi ■ ( 2 -12) 
The main result is that Qn(ji) can be expressed as the infinite sum, the so called 'cluster expansion' as: 

Q N (v) = Nl £ n MT? '? P . ( 2 - 13 ) 

z — ' - 1 A m, ! 

where 

m = (mi, . . . , mjv) . 

It must be stressed that eas. (l2.11[l - H2.13[l are purely combinatorial and they apply both for the short range interactions 
considered in ref.@ as well as the long range Newton potential. 

Since = O (4U for large N, and 5i,2,...,j contains at least a product of j factors fu, we see that 



Si 



2,...J 



, for N > j . 



Therefore, the cluster integrals ea. (|2.11[l take the form 

v y _1 



b j (v,N)={^} c 3 1 + 0(|) for /V»j. (2.14) 



N 



Here the coefficients Cj are positive numbers which only depend on the geometry of the box. As shown in fig. 1 the 
dominant terms for large N are cluster diagrams with a tree structure. In the large N limit, cluster diagrams with a 
loop structure as the last term in ea. (|2.12l) are subdominant. 

From eas. (l2~TU|) - (|2~T2) and lj2~Ti|) we find, 

ci = 1 , 



C2 



1 f 1 f l cPnd? 



r-2 



2 Jo Jo In - r 2 \ 



cs= 1 -[ 1 t f\f n fnl d3r i<- (2-15) 
2J 7 Jo \ri-r 2 \ \r 2 -f 3 \ 1 ' 



For a sphere of unit volume we find 



^sphere = 3 /4tt\ 7 = Q gg^ggj^j _ > = 51 /4tt\ 1 = ^32^3 _ _ _ ^ lg) 



For the cubic geometry chosen, it takes the value 2] 

c? be = 4 Al - a) /V - y) dy f . { \~ Z \ dZ = = 0.94116 . . . 

Jo Jo jo + v + ^ 

Furthermore, we find for the coefficient C4, 

c 4 = - I' I' I' t ^ d3r3 d3ri +- t t t t d3ri d3r2 d3r3 ^ (2 17) 



2 Jo Jo Jo Jo \ri -T2\ |f 2 -f 3 \ \r 3 -r 4 | 6 J J J Jo |n - r 2 | |f 2 - f 3 | |f 2 - f 4 
The combinatorial factors 5 and i take into account the symmetries of the respective cluster diagrams. 
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Inserting here the expansion in spherical harmonics, 



l m—+l 



where r> = max(r, r') and r< = min(r, r') , the angular integrals and then the radial integrals can be performed with 
the result, 

1 f 1 f 1 f 1 d 3 n d 3 r 2 d 3 r 3 d 3 r 4 _ 4 f b f b Uli^ &n _ 4tt 62 

JO JO JO JO l' r l ~ r 2| \''2 -T 3 | |T 3 -7-4| 

and 



o Jo 7 7 In ~r 2 | |r 2 -r 3 | |f 3 -r 4 | 7 Jo rf 2 r£ 3 r? 4 3 35 



Therefore, 



1 f 1 f 1 f 1 d 3 n d 3 r 2 d?r 3 rf 3 r 4 = f b f b jtf =1 rf d n = 4vr 188 

o Jo Jo Jo In - r 2 | |f 2 - r 3 | |f 2 - r 4 | J J r> 2 r> 3 r> 4 3 105 



sphere 47T 373 

° 4 " T 315 



(2.18) 



Notice that all these 'tree' cluster diagrams eas. J2.150 and i|2.17[l have a finite limit for zero cutoff i 



Divergent pieces for a — > are subdominant as for large N. The leading divergent contribution for a — > to 
Qn(v) to the nth. order in rj takes the formQ 

•r 1n ( n - 1) / ^ I =*- for " " 3 ' (2 . 19) 



n! TV™ 2 v ' y |fi -f 2 |™ 



log a for n = 3 . 



This gives for the physical quantities (see next section) contributions of the order 

Tj n 77 3 

— T a 3- " for n > 3 and — ^ logo for n = 3 . 

These contributions are clearly negligible in the N — > oo limit with fixed short-distance cutoff. 



III. THE MAYER EXPANSION FOR THE SELF-GRAVITATING GAS 

It is convenient to consider the generating function , 



X( V ,z) = J2 b MN)z^ , (3.1) 

3=1 



It can be shown that |6j 



J2 QjV ^ Z N = e X( V ,z) f (3 2) 

N=0 



where z is an auxiliary variable whose physical meaning (the fugacity) will appear in sec. IV. 
Hence, we can compute the coefficients Qn(v) from ea. (|3.2|l by contour integration, 

Q N (r)) r dz e x ^ 



N\ J 2m z^ 1 
We choose as contour a cercle of radius r, 

z = r e w , < 9 < 2tt 
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Integrating over 6 yields, 



Qn(v) 
M 



+7T 



(10 

2^ 



,I(i|,re'VN logi — iN i 



For large N, we can use eas. H2.14fl to express the bj(rf, N) in X(r), z) [see ea. (|3.1fl ] and we find, 



z rj^i N ( z r) s 



where 



Thus, 



3 = 1 



ivW JV2>i f +n d0_ w [i ff (r^%)_iogr-ie] 



AT! 



2vr 



(3.3) 



(3.4) 



(3.5) 



We can now apply the steepest descent method to this integral for large N since the integrand has the structure 

1 ( r e iB n\ 
*(r,0) = - ff (— -logr-ie, 

— (r,6)= l — 9 

(^v-J-r- (3 - 6) 

The saddle point, solution of ^(r,8) — = ^-(t*, 0) is thus found at 

M)saddle=(^0)> 
where t is A^-independent and is a solution of the equation, 

t g'(t t?) = 1 . (3.7) 

That is, t is a function of 77 denned by the constraint eq.J22J or more explicitly, 

+00 

j Cj (t rj) 1 = r) . 

3=1 

Choosing r = N t as integration path in ea. H3.5f) and expanding the integrand around 9 — yields, 



TV! ~ 



+ 00 



d6> 
27 



e JV *(JVt,0) 



l + O 



where, 



*(JVt,0) = -5(7? t„)-log[7V i„] 



a 2 * 



(Nt,0)=J2f Cj ti >0 



(3.8) 



where is a function of 77 defined by eq. (|3.7|) . 
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Using now Stirling's formula for the N\ factorial we find, 



N - ' •■■ 77 
Therefore, the free energy can be written as, 

F-F 



4 log'QwO?) = ~ g(v *v) - log*f, - 1 + O (jf 



(3.9) 



f 



vr =--g{vtr,) + logt n + l, 



where Fq stands for the free energy of the ideal gas, 

F = -NT log 

The pressure (at the surface) follows from the thermodynamic relation 

P = 

we find 



eV fmTY 
N \ 2^J 



dF 
dV 



(3.10) 



/'V)^^ = | + ^ and llogQ JV (,)=3 



NT 3 



3 Tj 



(3.11) 



The dimensionless ratio ■Stb was called f (rj) in refs.0, Q. We can express all physical quantities in terms of the 
function f(rj). We find from cas. (|3.9l) and l|3.11|l . 



3m = 2+ ihM and 



g(r) tn) + log/;,, + 1 = -3 



dx 



i - m 



That is, 



1 f v 
-logt v = f(r)) - 1 - dx 



1 - f{x) 



(3.12) 



(3.13) 



We can solve ea. <|3.7ll for t v in powers of r\ with the result, 

g(r) t„) =r]t n + c 2 (77 t n f + c 3 (77 Uf + O^ 4 ) , 

i„ = 1 - 2 c 2 V + [8 (c 2 ) 2 - 3 c 3 ] ?? 2 + 0(t7 3 ) . 
and from eq. (|3.12|l we find for f{rf) and £7(77 t v ), 

f{v) = 1 - I »7 + § [2 (c 2 ) 2 - c 3 ] t? 2 + ©(ry 3 ) , 

" da; l^M = _ C2 „ + [ 2 ( C2 )2 _ C3 1 ^ + 0{rj3) 



(3.14) 



3(77 t„) = 77 {1 - c 2 77 + 2 [2 (c 2 ) 2 - c 3 ] T7 2 } + 0(r? 4 ) 



(3.15) 



It must be noticed that the function / (77) (and hence all physical quantities) have the same expression whether we 
compute it from the mean field approach 0,0 or from the saddle point eas. (|3.5(l . I|3.9() and l|3.11|l . 



IV. THE GRAND CANONICAL ENSEMBLE FOR THE SELF-GRAVITATING GAS 



The definition of the variable z through ea. H3.2(l suggests that z is related to the fugacity. To be more precise, the 
grand partition function is defined in terms of the canonical partition function as, 



+ OC 



+00 



N=0 



N=0 



Qn{v) 



( mT 
V~2tT 



V z 



N 



(4.1) 



where we used ea. (|2.5|) . 

3 

Therefore, z must be multiplied by the ideal gas factor (^-) 2 V. 

At the saddle point we have after this renormalization 

I mT\ 2 If mT\ 2 
V 7. = N t„ and X I n.l 



, o J Vz = Nt v and X ^, ^— j V z j = logZ GC (v, z) ■ (4.2) 

Now, we know from ref.0] that the chemical potential /i takes the form 

' r & WM_3 [1 _ /(j7)] (43) 



T 

where /io stands for the chemical potential of the ideal gas, 



Ho = -T log 



V (niT 
N \~2^ 



(4.4) 



From ea. l|3.13|l and l|4.3|l we have, 



Therefore, from eas. l|4.2[) - (l4.5l) . 



^ = lo g t n . (4.5) 



t v = e t z , e T =z, (4.6) 

that is, we can identify z with the fugacity. 

We have found in ref. 2] that f(rf) obeys in the spherical case the first order differential equation of Abel's type, 

V R [3 f( V R ) - l]f{if) + [if {if) - 3 + V R ]f(v R ) = ■ (4.7) 
where the variable r/ R is defined as, 

4^ 1/3 



jf = i] = 1.61199. .. 7] . 

Using eas. (|4.7fl and (|3.13|l yields, 

log t„ = log , i.e. t n = e-i* f(r, R ) . (4.8) 

Ea. l|4.1(l provides the grand partition function. The above results showing that Zqc(j],z) is dominated by the 
canonical ensemble, together with eg. 14. 61) prove that the canonical and grand canonical ensembles are equivalent in 
their common region of validity as stated in refs.0,0]. 

In the thermodynamic limit, Zgc(Vi z ) IS given by eas. (|3.3|) and (|4.2|) as 

, ~ / x N fzn (mT\% \ N , , 

] og Z G c(Ti,z) = -g\^ J H=- fl (^t,). (4.9) 

Using here ea. l|3.12[) gives, 

logZ GC (v,z) = N[3f(T 1 )-2] , 
which exactly coincides with the expression found in ref. |2| for the grand partition function in the mean field approach. 



10 



V. CALCULATION OF THE CLUSTER COEFFICIENTS c 



We compute in this section the cluster coefficients c„ for the sphere by using a non-linear differential equation for 
the function g(u). 

We first find from ea. (|3.7f) 

/A \ 1/3 

ug' R (u) = 'n R with u = t n r] R and g R (u) = ( — J g(u) (5.1) 

Then, combining eas. (f4.7[) and If5.1f) yields the second order differential equation for g R (u), 

[u 2 g' R {u) + ug' R {u)} [2ug' R (u) + g R (u)] - 2 u 2 g R (u) - u g' R (u) + g R (u) = . (5.2) 

Alternatively, if we choose r) R = ug' R {u) as variable [see ea. (|5.1J) ]. we find from ea. (|4.7(l the first order non-linear 
differential equation 



[rf + 9n]^+V R [2v R + 9n-2]=0. 



That is, the invariance of ea. (|5.2|) under the rescaling of the variable u allows to reduce by one the order of the 
differential equation. 

Ea. H5.lfl has as regular solution around u = 0, 

g R (u) = J^c R u n with c%> here = c R , (5.3) 

and we can choose c R = 1. 

By inserting ea. l|5.3|l into ea. (|5.2|) the following nonlinear recurrence relation for the c R coefficients is obtained: 



c R 



" ( 2n +l)(„-l)^ 

We find, 



J2 c s ' c n- s s 2 (2n- 2s + 1) for n>2. (5.4) 



51 R _ 373 R _ 14911 R _ 2047 



r n _ i « _ _ -ft _ _ it _ _ - -f 

1 ~ ' 2 ~ 5 ' 3 ~ 70 ' 4 ~ 315 ' 5 " 6600 ' 6 ~ 429 ' ' ' ' 
in agreement with eqs. Q2.15[) . I|2.16|l and 12.1811 . 

By evaluating numerically the c R from the recurrence relation ea. H5.4|l . we find for large n: 



c R n ^ 0.309360.. n 5 



(u G c) n 



1 + 



(5.5) 



Here, uqc — 0.30034... is the radius of convergence of the series ea. (|5.3|l . Therefore, u — uqc is the nearest 
singularity of g R (u) in the u-plane. 

As we see from ea. l|4.9ll . g R (u) gives the grand canonical partition function as a function of u — t n r) R . Since 
u = t v r] R is proportional to the fugacity [see eq. I|4.6|l ] , uqc must be related to the critical point of the selfgravitating 
gas in the grand canonical ensemble. Indeed, using the critical value for rj R in the grand canonical ensemble 2], 

V R C = 0.79735 .. . and f( v R c ) 



3^c ' 
from ea. H4.8f) we obtain, 

«gc = Vgc tvga = Vac e~" Sc f(v§c) = \ = 0-30034 . 
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in perfect agreement with the value for uqc m ea. (|5.5() . The large order behaviour of the expansion coefficients 
eq. (|5.5|) corresponds to a (ugc ~ u ) 5 behaviour of the function gn(u). More precisely, for u — > uqc from ea. (|5.2|) we 
find, 



9R (u) U ^ GC 2(l-ry« c )+^ c 



UGC 



2 - vSc 1 



+ [(li- UGC ) 5 



Expanding this asymptotic behaviour in powers of u yields 

= 2(1-?7 GC )+77 GC . 



- 1 



E 



fe! u%t~, 



2 > u k 



This implies for the coefficients the following large order behaviour, 



7? n>l 



il 2-lcc »" 



7T (IGC)" 



1 + 01 - 

n 



which exactly coincides with ea. l|5.5|l since, 



1 /2 



0.309360 . 



In summary, the investigation here on the cluster expansion for the selfgravitating gas shows: 

• The selfgravitating gas admits a consistent thermodynamic limit N, V — > oo with -t- fixed. In this limit, 

extensive thermodynamic quantities like energy, free energy, entropy are proportional to N. That is, in this 
limit the energy, free energy, and entropy per particle are well defined and finite. 

• The cluster expansion and the mean field approach provide the same results in the thermodynamic limit 
AT, V ->■ oo with fixed. 

• The partition function needs a small short-distance cutoff a in order to be well defined, the thermodynamic 
limit has a finite limit for a — * 0. In other words, the contributions to the partition function that diverge for 
a — > are subdominant for N — > oo [see ea. (|2.19|) ]. That is, the N, V — oo limit of the cutoff model (with 

V 3 

fixed) has a finite limit for a — > 0. 



VI. THE STABILITY OF THE THERMODYNAMIC LIMIT N, V = L 3 -> oo WITH £ FIXED (FIXED 77) 

As shown in refs. 2] there are two phases in the selfgravitating gas: r/ < rj and 77 > 770 where 770 = 1.51024 . . . for 
spherical geometry and 770 ~ 1.515 for cubic geometry. For rj > 7/0 the selfgravitating gas collapses into a extremely 
dense phase with large and negative pressure. For 77 < 770 the selfgravitating gas is perfectly stable. The mean field 
applies in this gaseous phase and coincides with the expansion in powers of 77 discussed in the previous section. Such 
expansion converges within the gaseous phase. The variable 77 is related to the Jeans' length of the system as 

" = 3 (£) a ' L = Yl 

where 

[yf 1 _ n 

dj = \ in ' P = 17 

V m y/G m p V 

The gas collapses in the canonical ensemble for 77 > 770 which corresponds to L > dj. This corresponds to the Jeans' 
instability. 
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The relevance of the ratio " has been noticed on dimensional grounds[2j. However, the dimensionality argument 

alone cannot single out the crucial factor N in the variable 77. Notice that 77 contains the ratio and not 

y 3 v 

Therefore, in the thermodynamic limit 

V~N 3 and p=^~-L^0. 

As N, V — > 00, 77 is kept fixed in the same way as the temperature T. The energy, the product PV, the free energy, 
the entropy are expressed as N times functions of 77. The chemical potential, the specific heats, the compressibilities 
are just functions of rj. 

In a recent e-print it was stated that the thermodynamic functions for the selfgravitating gas diverge in the 
thermodynamic limit AT, V — ► 00 with fixed. 

We show here below that the statements made in ref. 7'j have crucial failures which invalidate the conclusions given 
in ref. |jj . 

Such statements in ref.Q are based in the inequality [eq.(30) in ref. 7]], 

1 r N v N 

Z C > W ^~ m/2 / n Y[d 3 n exp[-/35>,] > ^-expH^^J , 

v o 1=1 i<j ' i<j 

where 

/v 

1 f G 

(^j)vo = T7n II d3rk fa ' and fa = _ 1 ' 
v o J v N fe= i l r » r J'l 

and the last inequality follows from the property of the exponential function: (exp(y)) > exp((j/)). 

To write such inequality the author of ref. considers N selfgravitating particles in a volume V — R 3 with TV ~ R, 
that is, in the conditions of the above thermodynamic limit (i. e. N/V 1 ^ 3 fixed). 

Then, a portion of such volume of linear size Rq < R is considered: at this precise point the author assumes that 
./V r~j i?Q, that is, he assumes a extremely dense distribution of particles within the volume of size Rq. 

But for large N such distribution necessarily collapses since the gravitational gas avoids collapse only when 
N ~ Rq. In terms of the parameter 77, the assumption N ~ Rq implies 

770 - N/R ~ (i? ) 2 » 1 , 

which is deep in the collapsed phase |2j. That is, the assumption N ~ i?p necessarily implies that the gas collapses. 
The subsequent statements made in ref.0 are direct consequences of this assumption and are not valid. 



It is clearly true that the partition function sums over all configurations including collapsed situations. However, 
these collapsed configurations have a negligible weight for 77 < 770 ~ N/Rq. The argument of ref.Q applies only 
for 77 > 770 and only for such values of 77 collapsed states dominate the partition function. It must be noticed that 
in Nature, if collapsed configurations would always dominate selfgravitating systems, then, stars, galaxies and the 
interstellar medium would have collapsed since longtime. Fortunately, this argument in ref. 7] is not valid. 

As stressed in refs0, 0] there are two regimes for the selfgravitating gas in the canonical ensemble: 77 < 770 and 
77 > 770 with 770 = 1.51024 . . . for spherical geometry and 77c — 1.515 for cubic geometry. For 77 > 770 the selfgravitating 
gas do collapse into a extremely dense phase. It is to this collapsed phase and only to it that the eqs.(30)-(33) 
of ref. apply. 

The same comments applies to the convergence of the series expansions of the thermodynamic quantities discussed 
in ref. [eq.(34) of ref. above]. The series in powers of 77 only converge for q < rjc ■ Actually, the calculation of the 
radius of convergence of such series in ref. HQ provided us an independent check of the numerical value of 77c. 

In the same token, the thermodynamic limit proposed in ref. 8j leads to a cataclysmic collapse for the selfgravitating 
gas. It is proposed in ref. 8] to take N — > 00 with V ~ — > 0. 

It is obvious that in such limit the gas collapses since the volume vanishes. 
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(IV) f 12 f 23 f 34 = ■ ■ ■ ° 



(V) f 12 f 23 f 24 = 




FIG. 1: Diagrams contributing to the cluster expansion. Diagram (I) gives C2. Diagram (II) gives C3. Diagrams (III), (IV) and 
(V) contribute to a but (III) is subdominant for large N. 

More precisely, that means 77 ~ N's — > +00 deeply in the extremely dense phase. Furthermore, it is wrongly stated 
in sec. 4 of ref.[g] that the the entropy is not proportional to N as in refs.0,0. 
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